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ABSTRACT 

This  paper  presents  a  novel  technique  to  estimate  all  the  trajectory  parameters  of  fast  moving  targets 
using  a  single  SAR  sensor  without  increasing  the  PRF.  The  basic  reasoning  is  that,  although  the  returned 
echoes  may  be  aliased  in  the  azimuth  direction,  their  phase  and  amplitude  are  informative  with  respect 
to  the  moving  target  trajectory  parameters.  The  proposed  method  samples  the  data  in  the  spatial  domain 
and  extracts  the  data  along  the  signature  curve  that  depends  of  the  moving  target  trajectory  parameters. 
The  resulting  algorithm  is  a  highly  efficient  (from  the  computational  point  of  view)  ID  matched  filter. 
The  effectiveness  of  the  proposed  algorithm  is  illustrated  with  simulated  and  real  data. 

1  INTRODUCTION 

Unmanned  Air  Vehicle  (UAV)  systems  usually  benefit  from  low  weight  payloads.  UAVs  are  used  in 
the  battlefield  with  SAR  sensors  that  can  be  used  to  detect  moving  targets  and  estimate  their  velocity. 
The  problem  of  estimating  the  full  trajectory  parameters  of  moving  targets  using  SAR  data  typically 
requires  two  or  more  sensors.  Herein  we  present  a  novel  technique  to  estimate  the  full  velocity  vector 
of  moving  targets  that  is  light  from  the  computational  point  of  view  and  has  low  hardware  requirements 
since  it  uses  data  from  a  single  SAR  sensor.  The  requirement  of  a  single  SAR  sensor  contributes  to  the 
decrease  of  the  UAV  payload,  or,  alternatively,  enables  the  incorporation  of  other  sensors  thus  enhancing 
the  utility  of  the  overall  system. 

In  SAR  imaging  it  is  well  known  that  a  moving  target  induces  a  Doppler-shift  and  a  Doppler-spread  on 
the  returned  signal  in  the  slow-time  frequency  domain  [1].  Most  techniques  proposed  in  recent  literature 
to  focus  moving  targets  and  obtain  their  velocity  parameters  take  advantage  of  this  knowledge  [2],  [3]. 
The  azimuth  velocity  of  a  moving  target  is  responsible  for  the  spread  in  the  slow-time  frequency  domain, 
whereas  the  slant-range  velocity  is  responsible  for  the  Doppler  shift.  Given  a  PRF,  the  Doppler-shift, 
fD  =  2vx/X,  is  confined  to  fD  <  | PRF / 2 1 ,  where  vx  is  the  slant-range  velocity  of  the  moving  target  and 
A  is  the  carrier  wavelength.  If  the  induced  Doppler-shift  exceeds  PRF/2  (the  Nyquist  limit)  it  has  been 
mostly  accepted  that  the  true  moving  target  slant-range  velocity  cannot  be  uniquely  determined  using  a 
single  antenna  and  a  single  pulse  scheduling  [4] . 

In  [2]  we  have  already  presented  an  accurate  method  to  detect  and  estimate  all  the  moving  target 
parameters  using  a  single  SAR  sensor.  This  paper,  which  elaborates  on  ideas  presented  in  [5],  presents 
a  novel  technique  that  is  much  lighter,  from  the  computational  point  of  view,  to  estimate  the  same  pa¬ 
rameters.  Basically,  the  technique  corresponds  to  a  ID  matched  filter,  that  uses  the  signature  acquired  in 
the  spatial  domain  along  a  curve  that  depends  of  the  moving  target  trajectory  parameters.  By  using  both 
the  phase  and  the  amplitude  information  of  the  signature,  the  technique  is  able  to  retrieve  all  the  moving 
target  trajectory  parameters  using  a  single  SAR  sensor. 

We  start  by  writing  the  received  signal,  after  pulse  compression,  in  the  spatial  domain.  We  then  write 
explicitly  the  moving  target  signature  parameterized  by  the  trajectory  parameters.  By  assuming  that  the 
moving  target  signature  is  immersed  in  white  noise,  we  derive  a  maximum  likelihood  estimator  for  the 
moving  target  trajectory  parameters.  At  the  end  of  the  paper  we  present  results  using  simulated  and  real 
data  from  the  MSTAR  public  release  data  set. 
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Figure  1:  SAT?  slant-plane.  At  slow-time  u  0  t/re  radar  is  broadside  to  the  moving  target. 
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2  PROPOSED  APPROACH 

Let  us  consider  a  moving  target  with  slant-range  and  cross-range  coordinates  (x0.  y0),  when  the  plat¬ 
form  is  at  position  u  =  0,  and  velocity  ( vx ,  %)  =  (/ru#,  Aw#)  defined  in  the  slant-plane  ( x ,  y);  symbol  vR 
denotes  the  platform  speed  and  (p,  b)  is  the  target  relative  velocity  with  respect  to  the  radar.  The  moving 
target  coordinates  are  (x'  =  x0  —  pu0,  if  =  «o),  when  the  radar  platform  is  broadside  to  the  target,  as 
illustrated  in  Fig.  1. 

When  the  radar  is  positioned  at  coordinate  y  =  u,  the  corresponding  received  echo  from  a  moving 
target,  after  quadrature  demodulation  and  pulse  compression,  is  given  by 

u)  =  a{y0  -  vu)pc{x  -  r{u))fme~3‘2k°r{u\  (1) 


where  a(-)  is  the  two-way  antenna  radiation  pattern,  v  =  1  —  b.  pr ( x )  is  the  compressed  transmitted 
pulse  in  the  slant-range  direction,  Aq  =  27t/A  is  the  wavenumber  at  wavelength  A,  and  r ( u )  is  the 
distance  between  the  platform  and  the  moving  target,  given  by 

r{u)  =  \J{x o  -  pu)2  +  (y0  -  vu)2.  (2) 

It  is  assumed  that  the  dependency  of  a(-)  on  x  is  negligible. 

Let  us  dehne  v!  =  u  —  uq  such  that  u'  =  0  corresponds  to  the  position  of  the  radar  when  the  moving 
target  is  broadside  to  it.  Therefore,  defining  r'(u')  =  r(u'  +  u0),  we  obtain 


r'(u')  =  (x1  —  pu')2  +  {  —  vu')2. 


(3) 


Approximating  r'(u')  by  a  series  expansion  about  u! 
quadratic,  results 

r'(u')  ~  x'  —pu'  + 


=  0  and  retaining  only  the  terms  through  the 


(£+0  ,2 
2x' 


(4) 


valid  for  \u'\  x' .  If  the  range  migration  ip(u')  is  known,  then  we  can  form  a  vector  using  data  collected 

at  coordinates  s'(x,u')  =  s [x  =  x'  +  if(u' ),u'  +  u0\,  leading  to 


s'(x,  u') 


a(—uu')pc(0)fme 


j2k0r'(u') 


(5) 
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Since  pc(()  exhibits  high  resolution  about  (  =  0,  then  s'(x ,  u ')  becomes  clustered  about  x  =  x'  +  xp(u'), 
for  all  u'  in  the  support  of  a(-).  Therefore,  we  can  form  a  vector  containing  the  signature  samples  echoed 
by  the  moving  target. 

Let  us  assume  white  clutter  in  the  spatial  domain,  and  that  pc(0)  =  1.  Let  z(x,u')  =  s'(x,u')  + 
7i(x ,  u')  be  the  signal  plus  noise  returns.  Define  vectors 

Z  =  [Z-N,  .  .  .  ,Z0,  •  ■  ■  ,ZN]T 

S  =  [S-N,  •  ■  ■  5  So,  •  •  •  ,  Sn]T , 


where,  for  i  =  —N, . . . ,  N, 


Zi  =  z{x'  +  'lp(u'i),u'i) 

Si  =  a{  —  vu' i)e~^korl<'u'i\ 


with  u'i  =  u'si,  where  u's  is  the  sampling  space  in  the  cross-range  direction  chosen  such  that  the  clutter 
samples  exhibit  low  correlation. 

Notice  that  the  observation  vector  z  should  be  formed  with  data  sampled  along  the  curve  described 
by  xjj  and,  in  practical  situations,  data  is  available  only  at  discrete  values  of  ( x ,  u').  For  a  fixed  u'i,  the 
corresponding  slant-range  coordinate  is  computed  by  x'  +  ip  (u'i),  and,  in  general,  this  coordinate  does 
not  correspond  to  a  multiple  of  the  sampling  space  in  the  slant-range  direction.  Therefore,  a  resampling 
or  data  interpolation  should  be  done  in  the  slant-range  dimension.  Nevertheless,  as  illustrated  in  the 
Results  section,  a  simple  nearest  neighbor  interpolation  method  may  lead  to  satisfactory  results. 

Let  us  assume  that  the  moving  target  parameters  6  =  (//.  v.  x' .  y'.  i/0),  and  the  reflectivity  fm  are 
known.  In  this  case  only  the  noise  term  is  random.  Therefore,  the  density  of  vector  z  conditioned  to  6 
and  fm  is  Gaussian  given  by 

p(z|/m,  0)  -  N(nz,  cz),  (6) 


where  the  mean  is  nz  =  fms(6),  the  covariance  is  Cz  =  er^I,  and  the  clutter  cross-section  is  Notice 
that  Zi  =  Z;, (0) .  Consequently,  the  density  in  (6)  is  correct  only  for  the  true  6.  However,  for  0  ~  Q().  the 
range  migration  is  known  and  the  data  used  to  form  vector  s  is  taken  along  the  correct  curve  thus  leading 
to  z(G)  ~  z(6 o).  Thus,  we  still  use  (6)  to  compute  the  maximum  likelihood  (ML)  estimate  of  0.  After 
some  algebra,  we  obtain 


G 


ml  =  arg  max 

0 


£,^,•(9)1 

£,M0)|2  ' 


(7) 


We  note  that  ?/0  is  not  exactly  known.  However,  assuming  that  ?/ 0  is  a  multiple  of  us,  i.e.,  ?/0  =  usis, 
then  an  error  on  ?/0  means  an  integer  shift,  i0,  on  the  sequence  zt.  For  this  reason,  we  replace  (7)  by  the 
correlation 


(yU,  F,  =  arg  max 


max 

*0 


X,  (9)1) 

£,M9)|2  /• 


(8) 


Parameters  ?/0,  x0,  and  y0  are  then  given  as  follows: 


iio  =  u0  +  i0us ,  (9) 

x0  =  x'  +  fiu0 ,  (10) 

yQ  =  (z>-l  )u0.  (11) 
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The  proposed  methodology  can  be  summarized  as  follows: 

1.  Obtain  an  image  focused  with  static  ground  parameters.  The  ground  will  be  focused  and  the  moving 
targets  will  appear  defocused  and  at  wrong  positions. 

2.  Detect  the  moving  targets  and  estimate  their  approximate  positions  using  methods  such  as  those  pro¬ 
posed  in  [6]. 

3.  Isolate  the  defocused  moving  targets  and  re-synthesize  their  signatures  back  to  the  ( x ,  u)  domain  as 
proposed  in  [7].  This  procedure  increases  the  signal  to  noise  ratio,  although  it  introduces  some  correla¬ 
tion  due  to  the  spatial  windowing. 

4.  For  each  re- synthesized  signature  proceed  as  follows: 


•  Using  the  re- synthesized  signature  mass  center,  determine  a  rough  estimate  of  the  target  coordi¬ 
nates  (x1,  u0); 

•  Change  the  cross-range  coordinate  reference  according  to  u'  =  u  —  u0; 

•  For  all  (fi,  v)  pairs  of  interest,  form  the  observation  vector  z  using  data  extracted  along  the  curve 
given  by  (4)  and  compute  the  correlation  (8).  Estimate  as  (/i,  v)  the  velocity  pair  for  which  the 
maximum  absolute  value  occurs  and  store  its  coordinate  i 0. 

•  Estimate  (x0,  Vo)  using  equations  (9),  (10),  and  (11). 


Fig.  2  illustrates  the  overall  scheme  just  described,  in  a  simplified  scenario  containing  a  single  moving 
target.  On  the  top-right,  in  Fig.  2a),  the  MTI  contour  is  obtained  by  imaging  the  echoed  signal  after 
filtering  the  static  ground  returns  in  the  2D  frequency  domain.  The  coordinates  of  the  maximum  absolute 
value  are  estimated  as  the  position  (Xc,  Yc)  of  the  defocused  moving  object.  Using  this  coordinate  pair 
the  moving  object  is  digitally  spotlighted  in  the  spatial  domain  and  its  signature  is  resynthesized  back  to 
the  spatial  domain  as  illustrated  in  Fig.  2b).  After  estimating  the  signature  mass  center,  the  likelihood 
function  of  the  moving  target  velocity  is  computed,  as  shown  in  Fig.  2c).  Using  the  resynthesized 
signature  and  the  previously  estimated  moving  target  trajectory  parameters,  the  moving  object  is  then 
focused  and  repositioned  in  the  target  area  as  presented  in  Fig.  2d). 


3  RESULTS 

The  scheme  developed  in  the  previous  section  to  estimate  the  moving  targets  parameters  is  now 
evaluated  using  synthetic  and  real  data. 

A.  Synthetic  Data 

We  use  herein  three  types  of  simulated  targets:  point-like  targets,  extended  homogeneous  targets,  and  ex¬ 
tended  targets  with  predominant  scatterers.  All  the  targets  have  initial  coordinates  (x0,y o)  =  (64,  209)m, 
velocity  (vx.  vy  )  =  (—7.959,  8)m/s.  Notice  that  vx  induces  a  Doppler-shift  3  times  greater  than  the 
Nyquist  limit.  The  mission  parameters  used  in  the  simulations  are  summarized  in  Table  1. 

Table  2  presents  the  trajectory  parameter  estimation  results  for  the  point-like  target  as  function  of  the 
SCR  for  64  Monte-Carlo  runs.  These  results  show  that  if  the  SCR  is  higher  than  lOdB,  and  assuming 
that  the  targets  are  detectable,  then  the  algorithm  allows  the  velocity  and  initial  coordinates  estimation 
with  accuracy  enough  for  most  applications. 
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Target  region  in  the  2D  frequency  domain 
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Figure  2:  Schematic  summary  of  the  proposed  method. 
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Table  1:  Mission  parameters  used  in  the  simulation. 


Parameter 

Value 

Carrier  frequency 
Chirp  bandwidth 
Altitude 
Velocity 

Look  angle 

5  GHz 
100  MHz 

12  km 

637  km/h 
20° 

Table  2:  Trajectory  parameters  estimation  as  function  of  the  SCR. 


SCR  [dB] 

<Jvx  [m/s] 

crVy  [m/s] 

\yo-yo  [m] 

5 

0.18 

0.4 

65 

10 

0.04 

0.028 

3 

15 

0.01 

0.02 

1.4 

20 

0.004 

0.0013 

0.45 

Table  3  summarizes  the  estimation  results  for  the  three  considered  targets.  The  signal-to-clutter  ratio 
(SCR)  was  set  to  20dB.  As  can  be  seen,  if  the  target  is  point-like,  or,  if  it  has  some  predominant  scatterers, 
the  proposed  scheme  performs  quite  well  and  enables  the  velocity  estimation,  focusing  and  repositioning 
with  accuracy  high  enough  for  most  applications. 


B.  Real  Data 

We  now  present  results  using  a  BTR-60  vehicle  with  simulated  movement  and  a  background  scene  from 
Hunstville- Alabama,  both  taken  from  the  MSTAR  data.  The  mission  parameters  are  presented  in  Table 
4.  The  SCR  is  approximately  20  dB.  The  velocity  of  the  BTR-60  is  (vx,  vy)  =  (16.58,  2)  m/s  (vx  induces 
a  Doppler-shift  12  times  greater  than  the  Nyquist  limit). 

Fig.  3  presents  the  target  area  focused  with  static  ground  parameters.  As  expected,  the  vehicle  appears 
misplaced,  blurred  and  defocused.  If  correctly  processed,  the  BTR-60  image  should  appear  positioned  at 
coordinates  (x0,  yf)  =  (90, 124)  m.  After  detection,  we  digitally  spotlighted  the  BTR-60  signature  and 
resynthesized  it. 

Fig.  4  displays  the  resynthesized  signature  and  the  estimated  (x',u0)  =  (87, 123)m.  Using  the 
estimated  coordinates,  we  then  proceeded  with  the  algorithm  to  estimate  the  remaining  parameters.  The 
estimated  velocity  and  initial  coordinates  of  the  BTR-60  are,  respectively,  (vx,vy)  =  (16.53,  6.2)  m/s 
and  ( x0 ,  y0 )  =  (87.6, 129)  m. 

Using  the  estimated  velocity  and  initial  coordinates  we  focused  and  repositioned  the  BTR-60  as  shown 
in  Fig.  5. 

Table  3:  Simulation  results  (SCR=20  dB). 


Target  type 

vx  error 

vy  error 

X()  error 

r/o  error 

point-like 

0,005m/s 

0,012m/s 

0,02 

0,45m 

7  x  12m  extended 

0,3m/s 

4,2m/s 

0,18 

28,3m 

7x  12m  extended 
(one  predominant) 

0,041  m/s 

0,29m/s 

0,07 

1,6m 
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Table  4:  Mission  parameters  (MSTAR  data). 


Parameter 

Value 

Carrier  frequency 
Chirp  bandwidth 
Altitude 
Velocity 

Look  angle 
Oversampling  factor 

9.6  GHz 
250  MHz 

12  km 

637  km/h 

15° 

2 

BTR-60  unfocused 


64  96 

Slant-range  [m] 


128 


Figure  3:  Target  area  focused  using  the  wavefront  reconstruction  algorithm  with  static  ground  parame¬ 
ters. 
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BTR-60  resynthetized  signature 


Figure  4:  Resynthesized  signature. 
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BTR-60  focused  and  repositioned 


32  64  96  128 

Slant-range  [m] 


Figure  5:  Focused  and  repositioned  BTR-60. 
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Table  5:  Estimation  results  for  BTR-60  (SCR=  20  dB). 
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Target  type 

vx  error 

vy  error 

xq  error 

Vo  error 

BTR-60 

0,005m/s 

4,2m/s 

2,4m 

5  m 

Notice  that  since  we  are  using  a  single  SAR  sensor  and  because  the  moving  target  spectra  (although 
6  times  folded)  is  completely  superpositioned  on  that  from  the  clutter,  we  cannot  remove  the  defocused 
BTR-60  from  the  original  image.  To  achieve  this  goal  the  moving  target  spectrum  should  not  overlap 
completely  the  clutter  spectrum.  Another  possibility  would  be  to  have  data  from  more  than  one  SAR 
sensor  [8]. 

The  obtained  results  for  the  BTR-60  under  study  are  summarized  on  Table  5.  As  can  be  seen,  the 
proposed  methodology  yields  velocity  and  position  estimates  accurate  enough  for  focusing  and  reposi¬ 
tioning  the  vehicle. 

From  the  results  presented  in  this  section,  we  conclude  that  the  suggested  strategy  works  well  for 
point-like  targets  or  extended  targets  exhibiting  some  predominant  targets.  When  this  is  not  the  case, 
the  algorithm  still  gives  good  results  for  the  slant-range  velocity.  It  produces,  however,  large  errors  on 
the  cross-range  velocity  estimation  and  on  the  estimated  cross-range  initial  position.  Nevertheless,  the 
presented  methodology  can  be  used  to  estimate  the  velocity  and  initial  coordinates  of  most  man-made 
targets,  since  they  typically  exhibit  predominant  scatterers. 


4  CONCLUSIONS 

In  this  paper  we  presented  a  spatial  based  methodology  for  the  estimation  of  all  the  moving  target 
parameters  using  a  single  sensor.  The  main  algorithm  samples  the  spatial  domain,  compressed  in  the 
slant-range  direction  and  uncompressed  in  the  cross-range  direction,  to  acquire  data  along  the  signature 
curve  that  depends  of  the  moving  target  kinematics.  To  achieve  efficiency  and  simplicity  we  derived  the 
ML  estimator  of  the  velocity  parameters  by  assuming  white  noise.  This  assumption  led  us  to  a  matched 
filter  type  solution.  The  method  was  successfully  tested  using  a  combination  of  simulated  and  real  data, 
with  all  targets  inducing  Doppler-shifts  beyond  the  Nyquist  bound  imposed  by  the  mission  PRF. 
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